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COMPUTER PROGRAMS FOR LOCATING THE MOON 


In two previous issues of the EIMAC EME notes, methods were described 
for determining antenna aiming data to "hit" the moon. In Bulletin AS-49, 
there was a reprint of an article by Don Lund, WA0IQN, an article by 
Victor Michael, W3SDZ, and some letters to the Technical Editor of QST 
from Alan Goodacre, VE3BZS, all discussing antenna aiming. Also, in Bulletin 
AS-49-1 more information on aiming data was presented. Bulletin AS-49-1 
stressed the use of the Nautical Almanac and the Hydrographic Office publi- 
cation HO-214. It is possible to determine the moon's position by means of 
these two publications without extensive mathematics. However, as one be- 
comes more involved in the EME project, the habits of the moon become 
interesting. Also, more sophisticated methods of tracking the moon become 
attractive. 

For those persons having access to the hand-held scientific calculators, 
new techniques become available. For those amateurs having access to a 
computer, either punch card, or an active terminal, another interesting 
hobby of writing and modifying programs develops. 

In this issue of the EME notes, four programs are presented which may 
be of interest. One is in BASIC, two are in TYMSHARE SUPERBASIC and one in 
TYMSHARE SUPERFORTRAN. As is usually the case, these programs will probably 
need to be modified to run on the machine available to the individual 
amateur. 

It should be emphasized that a computer is not necessary to have suc- 
cessful moonbounce contacts. A computer is just another interesting tool to 
work with which can add to the fun. 

PROGRAM I 
LANGUAGE: BASIC 

This program will run on almost all machines accepting BASIC. The 
operator must put in data from the Nautical Almanac for the year, month, 
day and hours that the elevation and azimuth is desired. Also, the latitude 
and longitude of the geographic location in question must be entered. State- 
ment 20 allows the operator to put in the time interval, in minutes, between 
each calculation. 

The origin of this program is not known. 

PROGRAM 2 

LANGUGAGE: TYMSHARE SUPERBASIC 

The program is a modified version of Program One. The output format 
is slightly different. Also, limits can be put on the azimuth and elevation 
angles at the operator's discretion. The program requires latitude and 
longitude and time interval inputs. Also, Nautical Almanac data for the 
year, month, day and hours must be put into the data statements. Statement 
610 allows the operator to put the data in for identification. The computer 
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does not use this statement in the calculations. The computer will pause 
and let you type in the date. When inserting the date, do not use commas 
for separating the month, day and year. If you do, the computer will become 
confused. When the carriage return is depressed the computer will carry 
on. 


PROGRAM 3 

LANGUAGE: TYMSHARE SUPERBASIC 

This program was originally written in BASIC by Lance Collister, 
WA3GPL. However, the BASIC used on the General Electric terminal was 
sufficiently different from the TYMSHARE SUPERBASIC, that the program 
would not run. A few modifications were made to allow the TYMSHARE machine 
to run. What is presented here will not run on the GE computer. Perhaps 
those who attempt to use this program will have to make further changes 
for their individual cases. 

The program is very useful because no almanac data is required. The 
program allows the calculation of the moon position by comparing back to 
the position of the moon at the start of the Julian calendar. The operator 
need put in only the latitude, longitude, time interval, month, day, year 
and hours between which aiming data is required. If rising, or setting, 
moon calculations are required, the computer determines the time and will 
execute the printout requested. 

PROGRAM 4 

LANGUAGE: TYMSHARE SUPERFORTRAN 

This program was supplied by Louis Anciaux, WB6NMT. As received, the 
program was on tape and written in Fortran IV. Again, the program was modi- 
fied to run with TYMSHARE SUPERFORTRAN. This program requires Nautical 
Almanac data. However, the GHA and DECLINATION for 0000, 0600, 1200, 1800, 
and 2400 hours GMT is all that is required. The program will interpolate 
for all other times. Of course the latitude, longitude and time interval 
for the calculation must be typed in. 
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PROGRAM I 
LANGUAGE: BASIC 


10 DATA 37* 23- 5, 122, 10*5 
20 DATA S 

100 REM THIS PROGRAM CALCULATES MOON AIMING DATA FROM A SPECIFIED 
110 REM OTH AT SELECTED TIME INTERVALS DURING THE HOUR* 

120 REM 

130 REM DATA REQUIREMENTS ARE AS FOLLOWS: 

140 REM 

150 REM LINE 10 - LATITUDE DEGREES* LATITUDE MINUTES* 

150 REM LONGITUDE DEGREES* LONG I DUDE MINUTES 

170 REM 

180 REM LINE 20 - TIME INTERVAL FOR CALCULATIONS IN MINUTES 
19 0 REM 

200 REM LI NF 1 050- -ALMANAC DATA CONSISTING OF THE TIME* GHA DEGREES 
210 REM GHA MINUTES* DECLINATION DEGREES AND DECLINATION 

220 REM MINUTES 

230 RFM 

240 RFM DATA MUST BE GIVEN FOR BOTH THE BEGINNING AND 

250 REM END OF THE HOUR. THAT IS* IF AIMING DATA IS 

250 REM DESIRED FOR HOURS BEGINNING 09 00* 1 000* AND 1100 

270 REM YOU MUST SHOW DATA FOR 0900*1000*1100 AND 1200. 

280 REM 

290 REM TO INDICATE THE END OF THE HOURLY DATA INCLUDE A SINGLE 
300 RFM DATA LINE WITH A ONE* A TWO FOR FINISHED MONTH 

305 REM WITH THIS METHOD A MONTH AT A TIME CAN BE DONE 

310 REM 

320 LET R= 1 . 74533E-02 
330 LET PI = 180* h' 

340 LFT P3= 350*R 

350 READ L 3* L4 

350 LFT L 1 =( L3+L4/60) *R 

370 READ L 3* L4 

380 LET L2=(L3+L4/50)*R 

390 LET L3=SINCL1) 

AQO LFT L4=C0SCL1) 

405 READ J1 
405 LET J= 60/ J 1 
410 PRINT 
420 PRINT 
430 PRINT 

440 PRINT "MOON AIMING DATA SPECIALLY PREPARED FOR K6MYC" 

450 PRINT 
450 PRINT 

470 PRINT "TIME”* "AZIMUTH"* "ELEVATION” 

480 PRINT 

510 READ B 

530 READ B 1 * B2* B3* B4 

540 LET G1=(B1+B2/60>*R 

550 LET Dl=( B3+B4/60)*R 

550 READ B5 

575 IF B5= 1 THEN 450 
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PROGRAM I (cont'd) 


576 IF B5-2 THEN 2000 

577 IF B5= 3 THEN 48 0 
58 0 HEAD B1*B2,B3*B4 
590 LET G2=CB1+B2/60)*R 
600 LET D2=( B3+B4/60) *R 
610 LET G0=(G2-G1)/J 
620 IF G0> 0 THEN 640 
630 LET G0=( C-2+P3-GD/J 
640 LET D0=(D2-D1)/J 
650 LET G=G 1 

660 LET D=D 1 

670 FOR 1=0 TO 59 STEP J1 
680 LET D3=SIN(D) 

690 LET D4=C0SCD) 

700 LET H=L2-G 
710 IF H-P 1< 0 THEN 750 
720 IF H-P 1 = 0 THEN 770 
730 LET H= F3-H 
740 GOTO 770 

750 IF H+P 1 >= 0 THEN 770 
760 LET H= P3+H 
770 LET A= 0 
78 0 LET E= 0 

790 LET E3=L 3* D3+L4* D4* COSC H ) 
800 IF E3< 0 THEN 940 
810 LET E4=S0R( 1-E3t 2) 

820 LET A3=SIN(H)* D4/E4 
830 LET A4=SQR( 1-A3t2) 

840 LET A= ATN( A3/A4) 

850 LET E0=E3/E4-4/C 240*E4> 
860 IF EO <= 0 THEN 950 
870 LET E=ATNCEO)/R 
880 LET E=INT< 10*E+. 55/10 
890 IF L3*E3-D3>0 THEN 930 
900 IF A >= 0 THEN 940 
910 LET A= P3 + A 
920 GOTO 940 
930 LET A= P 1 - A 
940 LET A=INT( 10+A/R+.5)/ 10 
945 IF E=0 THEN 960 
950 PRINT B+ 1 ,» A* E 
960 LET G= G+GO 
970 LET D= D+ DO 
980 NEXT I 
990 LET G 1 = G2 
1000 LET D 1 =D2 
1010 LET B=B5 
1020 GOTO 560 
2000 END 
> 
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PROGRAM 2 

LANGUAGE: TYMSHARE SUPERBASIC 


100 DATA 37* 34, 122, 18 
110 DATA 15 

120 REM THIS PROGRAM CALCULATES MOON AIMING DATA FROM A SPECIFIED 
130 PEM QTH AT SELECTED TIME INTERVALS DURING THE HOUR. 

140 REM 

150 REM DATA REQUIREMENTS ARE AS FOLLOWS: 

160 REM 

170 REM LINE 100 -LATITUDE DEGREES, LATITUDE MINUTES, 

180 REM LONGITUDE DEGREES, LONGITUDE MINUTES 

190 RFM 

200 REM LINF 110 -TIME INTERVAL FOR CALCULATIONS IN MINUTES 
2 1 0 REM 

220 REM LINE 1 270- -ALMANAC DATA CONSISTING OF THE TIME, GHA DEGREES 
230 REM GHA MINUTES, DECLINATION DEGREES AND DECLINATION 

240 PEM MINUTES 

250 REM 

260 REM DATA MUST BE GIVEN FOR BOTH THE BEGINNING AND 

270 REM END OF THE HOUR. THAT IS, IF AIMING DATA IS 

280 PEM DESIRED FOR HOURS BEGINNING 09 00, 1 000, AND 1100 

290 REM YOU MUST SHOW DATA FOR 0900,1000,1100 AND 1200. 

300 REM 

310 REM TO INDICATE THE END OF THE HOURLY DATA INCLUDE A SINGLE 
320 PEM DATA LINE WITH A ONE, A TWO FOR FINISHED MONTH 

330 REM WITH THIS METHOD A MONTH AT A TIME CAN EE DONE 

340 REM 

350 LET R= 1 . 74533E-02 

360 LET P 1 = 1 8 0*R 

370 LET P3= 360*R 

38 0 READ L 3, L4 

390 LET L1 = CL3 + L4/60)*R 

400 READ L 3, L4 

410 LET L2=(L3+L4/60)*R 

420 LET L3=SINCL1> 

430 LET L4= COS C L 1 ) 

440 READ J1 
450 LET J=60/U1 
460 PRINT 
470 PRINT 
480 PRINT 

490 PRINT "MOON AIMING DATA FOR W6P0." 
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PROGRAM 2 (cont'd) 


500 

510 

520 

530 

540 

550 

560 

570 

58 0 

59 0 
600 
610 
620 
630 
640 
650 
660 
670 
68 0 
69 0 
700 
710 
720 
730 
740 
750 
760 
770 

78 0 

79 0 
800 
8 10 
820 
8 30 
840 
850 
860 
870 
880 
890 


PRINT 

PRINT "MINIMUM AZIMUTH": 

INPUT Z1 

PRINT "MAXIMUM AZIMUTH": 

INPUT Z2 

PRINT "MINIMUM ELEVATION": 

INPUT M 

PRINT "MAXIMUM ELEVATION": 

INPUT N 

PRINT 

PRINT 

PRINT "DATE": 

INPUT X 

PRINT 

PRINT 

PRINT "GMT" * "AZIMUTH'S "ELEVATION" 

PRINT 

READ B 

READ B1,B2jB3jB4 
LET G1=(B1+B2/60)*R 
LET D 1 = < B3+B4/60 ) * R 
READ B5 

IE B5= 1 THEN 590 
IE B5= 2 THEN 2000 
IF B5= 3 THEN 660 
RFAD B 1 j B2> B3* E4 
LET G2=( B1+B2/60)*R 
LET D2= ( P3+B4/60) *R 
LET G0=CG2-G1)/J 
IF G0>0 THEN 810 
LET G0=( G2 + P3-G1 ) /J 
LET D0 = ( D2-D 1 ) / J 
LET G= G 1 
LET D=D 1 

FOR 1=0 TO 59 STEP J1 
LET D3=SIN(D) 

LET D4= COS CD) 

LET H=L2- G 

IF H- P 1 < 0 THEN 920 

IE H-P 1 = 0 THEN 940 
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PROGRAM 2 (cont'd) 


900 LET H= P3-H 
910 GOTO 940 

920 IF H + Pl >= 0 THEN 940 
930 LET H=F3+H 
940 LET A= 0 
9 50 LFT E= 0 

950 LFT E3=L3* D3+L4* D4*C0 SC H ) 

970 I E E3<0 THEN 1110 
980 LFT E4=SQRC 1-E3t 2) 

990 LET A 3= SI N C H > * D4/E4 
1000 LET A4=SQR< 1-A3t2) 

1010 LET A=ATN( A3/A4) 

1020 LET E0=E3/E4-4/( 240*E4) 

1030 IF EO <= 0 THEN 1 170 

1040 LET E=ATN(EO)/R 

1050 LET E= I NTC 1 0*E+ • 5 ) / 1 0 

1060 IE L3*E3-D3>0 THEN 1100 

1070 IE A >= 0 THEN 1 1 10 

1080 LET A=F3+ A 

1090 GOTO 1110 

1100 LET A=P 1- A 

1110 LET A= I NTC 1 0* A/R+ • 5 ) / 1 0 

1120 IE E=0 THEN 1200 

1130 IE A>Z 1 THEN 1140 ELSE 1200 

1140 IE A<Z 2 THEN 1150 ELSE 1200 

1150 IF E>M THEN 1160 ELSE 1200 

1160 IF E<N THEN 1170 ELSE 1200 

1170 PRINT IN FORK "4D":(B+I) 

1180 PRINT IN IMAGE " %%%.% %%.%"SAjE 

1 190 IF 1 = 45 THEN PRINT 

1200 LFT G= G+GO 

1210 LET D=D+DO 

1220 NFXT I 

1230 LET G 1 = G2 

1240 LET D 1 =D2 

1250 LET B=B5 

1260 GOTO 710 

2000 END 
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PROGRAM 3 

LANGUAGE: TYMSHARE BUPERBASIC 

100! THIS PROGRAM IS DESIGNED TO CALCULATE THE AZIMUTH AND 

1 10 {ELEVATION OF THE MOON. 

120 ! 

130! THE REQUIRED INPUT IS THE GMT MONTH* DAY* AND YEAR FOR 

140! WHICH THE MOON'S AZIMUTH AND ELEVATION ARE DESIRED. IN 
150! ADDITION* IF PRINTOUT IS DESIRED NOT ONLY WHEN THE MOON IS NEAR 
160! THE HORIZON* THE TIME INTERVAL OVER WHICH THE POSITION IS DESIRED 
170! MUST BE ENTERED. 

180! 

190! THE COMPLETE INPUT FORMAT FOR THE ABOVE INFORMATION IS: 

200 ! 

210! MM* DD* YYYY* TTTT* TTTT 

220 ! 

230! IF YOU ELECT TO HAVE PRINTOUT OCCUR ONLY WHEN THE MOON 

240 ! NEAR THE HORIZON* YOU WILL BE ASKED TO INPUT THE MAXIMUM ELEVATION 
250! FOR WHICH YOU WANT PRINTOUT- IN THIS CASE* YOU WILL NOT BE ASKED 
260! TO SUPPLY TIMES FOR THE TIME INTERVALS. 

270! 

280! PRINTOUT IS SUSPENDED WHENEVER THE ELEVATION OF THE MOON IS 

29 0 ! NEGATIVE • 

300! 

310! YOU CAN ASK FOR DATA FOR UP TO 25 DAYS AT ONCE: SIMPLY 

320 ! TYPE IN THE INFORMATION FOR ONE DAY CAS IN THE REQUIRED FORMAT 
330! FOR YOUR PARTICULAR NEEDS) EACH TIME THE COMPUTER PRINTS A ? • 

340 '.WHEN YOU HAVE ENTERED ALL THE DATA YOU WISH* FOLLOW THE LAST ? BY AN 
350! INPUT OF ZEROES SEPARATED BY COMMAS- JUST AS IN THE FORMAT YOU WERE 
360JUSING TO ENTER THE DATA. 

370! 

380! ALL DATES AND TIMES USED IN THIS PROGRAM ARE IN GREENWICH MEAN 

39 0 ! TI ME • IN ADDITION* ALL TIMES ARE IN THE 0000 TO 2400 HOUR SYSTEM. 
400! 

410! THE BASIC USED IN THIS PROGRAM HAS BEEN KEPT AS ELEMENTARY AS 

420! POSSIBLE TO FACILITATE ITS USE ON OTHER TIME SHARING SYSTEMS. I 
430IWISH TO ACKNOWLEDGE THE GENEROUS ASSISTANCE OF THOMAS AKE OF THE 
440 ! WARNER- S WAS EY OBSERVATORY* WITHOUT WHOM THE VITAL EQUATIONS FOR 
450! COMPUTING THE MOON'S POSITION WOULD NOT HAVE BEEN AVAILABLE. 


460! 


470! 

LANCE COLLI STER 

48 0! 

WA 1 JXN/WA3GPL 

49 0! 

CLEVELAND* OHIO 

500! 

MARCH* 1974 

510! 



520 ! ******************************************************************** 

530 DIM F( 25)*V(25)*Y( 25)* Q( 25) * SC 25) 

540 DEF FNAC X ) = I NTC X*D5* 1 0+ • 5 ) / 1 0 

550 DEF FNBCX)=CX-INTCX) >*P5 

560 LET P5=2. 0000000000*3. 1415926535 

570 LET D5= 360. 000000000 0/P5 ! CONVERSION TO DEGREES 
580 LET R5=P5/360. 0000000000 ! CONVERSION TO RADIANS 

583 STRING Z8 

584 PRINT 

585 PRINT "WHAT ARE THE CALL LETTERS OF THE STATION": 

586 INPUT Z8 

590 PRINT "WHAT IS THE LATITUDE IN DEGREES* MI NUTES" : 
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PROGRAM 3 (cont'd) 


600 INPUT L 5 * U 5 

610 PRINT "WHAT IS THE LONGITUDE IN DEGREES* MI NOTES": 

620 INPUT L 6 * U6 

630 LET L5=CL5+U5/60)*R5 

640 LET L6= ( L6+U6/60) *R5 

650 PRINT "WHAT IS THE DESIRED PRINTING INCREMENT IN MINUTES": 

660 INPUT I 

670 PRINT "DO YOU ONLY WANT PRINTOUT WHEN THE MOON" 

680 PRINT "IS NEAR THE HORIZON": 

69 0 INPUT B3> 

700 IF BS="YES" THEN 730 
710 LET 16=100 
720 GOTO 800 

730 PR I NT"BELOW WHAT ELEVATION IN DEGREES DO YOU WANT PRINTOUT TO OCCUR? 
740 INPUT 16 

750 PRINT "WHAT ARE THE GMT MONTH* DAY * YEAR DESIRED": 

760 FOR N= 1 TO 25 

770 INPUT FCN)*V<N) *Y<N) 

780 IF F( N ) =0 THEN 860 
78 5 NEXT N 
790 GOTO 760 

800 PRINT "WHAT ARE THE GMT MONTH* DAY * YEAR* TIME I NTERVALC BEGI NNI NG* " 

810 PRINT "ENDING) DESIRED": 

820 FOR N= 1 TO 25 

830 INPUT F<N) *VCN) *Y(N>* QCN)* SCN) 

840 IF F( N ) =0 THEN 860 

845 NEXT N 

850 GOTO 820 

860 LET N 5=N- 1 

870 FOR N= 1 TO N5 

880 IF B$="YES" THEN 900 

890 GOTO 930 

900 LET E 1=2400 

910 LET B=0 

920 GOTO 950 

930 LET E 1= S( N ) 

940 LET B=Q( N ) 

950 LET M=F(N) 

960 LET D=VCN) 

970 LET Y=Y ( N ) 

980 LET Y 1=Y- ( INT(Y/100)*100) 

990 PRINT 
1000 PRINT 

1010 PRINT "POSITION OF MOON ON " : M : "/" : D: "✓ " : Y 1 : "GMT" 

1020 PRINT 

1030 PRINT "GMT"* "AZ"* "EL"* "GHA"* "DEC" 

1040 PRINT " "*" — "*" — "*" ", " " 

1050 PRINT 
1060 LET I 1 = 2 

1070 !HERE BEGINS CALCULATION OF JULIAN DATE 
1080 IF M>= 3 THEN 1160 
1090 IF INTC CY- 1853)/4)< 1 1 THEN 1120 
1100 LET C 1 = - 1 
1110 GOTO 1130 
> 
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PROGRAM 3 (cont'd) 


1120 
1130 
1 140 
1150 
1160 
1170 
1 180 
1 190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
128 0 
129 0 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 

148 0 

149 0 
1500 
1510 
1520 
1530 
1540 
1550 


LET C 1=0 

LET J 1 = 365* ( Y- 18 53) +D+ 30* C M+9 ) +INTC C M+ 10) /2> 

LET J2=INTC (Y- 18 53) /4) + 1+C 1 
GOTO 1270 

IF INTC C Y- 1852)/4)< 1 1 THEN 1190 

LET C 1 =- 1 

GOTO 1200 

LET C 1=0 

IF M=9 THEN 1240 

IF M= 1 1 THEN 1240 

LET C2=0 

GOTO 1250 

LET C2= 1 

LET J 1=365* CY- 18 52) + D+ 30* ( M- 3) + 1 NT( CM-2)/2) 

LET J2=INTC (Y- 1852) /4) +C 1+C2 

LET J=U1+J2 !( JULIAN DATE- 239 7547 • 5) FOR 0 HOURS GMT 
LET T 1= J- 1 7472 • 5 

LET D9 = ( B-INTC B/ 1 00) * 1 00) + 1 NTC B/ 100) *60 

LET D6= CE1-INTCE1/100)* 100) + INTC El/ 100)*60 

LET D7=D9-D6 

LET B8=D7-I 

IF D7> 0 THEN 1350 

GOTO 138 0 

IF D8>=0 THEN 2220 
LET B=E 1 

! CALCULATION OF LATITUDE AND LONGITUDE OF MOON 
LET T= C B- 1 NTC B/ 100)* 100)/ 1440+1 NTC B/ 100)/ 24 
LET T 5 = T 1 + T 

LET K 1=FNBC . 75 1 2 1 3+ . 03660 1 1 02* T5 ) 

LET K2=FNBC . 8 225 1 3+ . 03629 1 6457*T5) 

LET K3=FNBC . 99 5766+ • 002737778 52*T5) 

LET K4=FNBC . 9 7427 1 + . 0338 63 19 22*T5 ) 

LET K5=FNBC . 03 1 2525+ . 0 36748 1957*T5) 

LET L8=K 1+. 658*R5*SINC 2*K4) +6 . 289*R5*SI NC K2) 

LET L8=L8- 1 • 274*R5* SI N C K2- 2*K4 ) - • 1 8 6*R5* SI NC K3) 

LET L8=L8+.214*R5*SINC 2*K2)-. 1 14*R5*SINC 2*K5) 

LET L8 =L8- • 0 59*R5* SI NC 2*K2- 2*K4 ) - . 057* R5* SI NC K2+K3- 2*K4 ) 

LET K6=K5+ . 6 59 3*R5*SINC 2*K4 ) + 6 • 2303* R5* SI N C K2) - 1 . 272*R5*SI NC K2- 2*K4 
LET L 7= 5 • 144*R5*SINCK6) - • 146*R5*SI NC K5-2+K4) 

! CALCTJLATI ON OF RIGHT ASCENSION CA=R1) AND DECLINATION CDl) 

LET D1=C0SCL7)*SINCL8)*. 3978 2 1+SINC L7) * -9 17463 
LET D1=ATNC Dl/C SORC 1-Dlt 2) ) ) 

LET A2=C0SCL7)*C0SCL8)/C0SCD1) 

LET A 1 = C COSC L7)*SI NC L8 )* -9 1746 3-SINCL7)** 397821) /CO SC Dl) 


1560 LET A= ATNC A 1 / A2) 

1570 GOSUB 1860 

1580 LET R 1 =A 

1590 LET Ll=. 065709822+T1 

1600 LET L=T*24*1. 002738 + 6. 646055+ C L 1 - 1 NTC L 1/24) * 24) 

1610 LET L=CL-INTCL/24)*24) 

1620! CALCULATION OF GREENWICH HOUR ANGLE* G* FROM LOCAL SIDEREAL TIME 
1630 LET G=CL/24)*P5-R1 
1640 IF G<P5 THEN 1670 
1650 G=G-P5 


> 
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PROGRAM 3 (cont'd) 


1660 GOTO 1710 
1670 IF G< 0 THEN 169 0 
1680 GOTO 1710 
1690 G=G+P5 

1700! CALCULATION OF YOUR LOCAL HOUR ANGLE,H, FROM GHA 
1710 LET H=L6-G 

1720! CALCULATION OF ELEVATI ON* E* OF OBJECT 

1730 LET E3=C0S<L5>*C0SCH)*C0S(D1)+SINCD1>*SINCL5> !SIN OF ELEVATION 
1740 LET E2=SQRC 1-CE3*E3) ) 

1750 LET E=ATN( E3/E2) 

1760 IF E<0 THEN 2170 
1770 IF E>I 6*R5 THEN 2170 
1780! CALCULATION OF AZIMUTH, A, OF OBJECT 
1790 LET A2=SIN( D1)/CC0S(L5)*C0S( E) ) 

1800 LET A2=A2-C SINCL5)/C0SCL5> )*C SINC E)/COS( E> ) 

1810 LET A1=SINCL5)*SIN(D1)+C0S(L5)*C0SCD1)*C0SCH) 

1820 LET A 1 = ( SINCH > *COS( D 1 ) ) /( SQRC 1 -A 1 1 2) ) 

1830 LET A=ATN< A1/A2) ! AZ I MUTH= ARCTANC SIN OF A Z/ COS OF AZ ) 

1840 GOSUB 1860 
1850 GOTO 2020 

1860! REMOVAL OF AMBIGUITIES INCURRED WITH ATN FUNCTION 

1870 IF A=0 THEN 1890 

1880 GOTO 1930 

1890 IF A2<0 THEN 1910 

1900 GOTO 2010 

19 10 LET A= P5/2 

1920 GOTO 2010 

1930 IF A>0 THEN 1990 

1940 IF A2<0 THEN 1970 

1950 LET A=P5+A 

1960 GOTO 2010 

1970 LET A=P5+C A-P5/2) 

1980 GOTO 2010 
1990 IF A2>=0 THEN 2010 
2000 LET A=A+ P5/2 
2010 RETURN 

2020 IF ( T-I 1 )>( 2*1 )/ 1440 THEN 2040 
2030 GOTO 2145 
2040 PRINT 

2145 PRINT IN FORM "4D" : ( I NTC B+ • 5 ) ) 

2150 PRINT IN IMAGE " %%%.% %%.% 

**%•* %%%•%": FNAC A) , FNAC E ) * FNAC G) , FNAC D 1 ) 

2160 LET I 1=T 
2170 LET B=B+I 

2180 LET Z = ( B“I NTC B/ 100)* 100)-60 

2190 IF Z<0 THEN 1290 

2200 LET B= INT( B/ 100)*100+100+Z 

2210 GOTO 1290 

2220 NEXT N 

2230 LET N=0 

2240 PRINT 

2250 PRINT 

2260 PRINT 

2270 PRI NT*'DO YOU WANT MORE INFORMATION": 

2280 INPUT D$ 

2290 IF D35="YES H THEN 540 
2300 END 


- 11 - 



PROGRAM 4 

LANGUAGE: TYMSHARE SUPERFORTRAN 


100 


101 


102 


103 

100 

104 

101 

105 

102 

106 

103 

107 

104 

108 

105 

109 

106 

110 

107 

1 1 1 

108 

112 

109 

113 

1 10 

1 1A 

1 1 1 

115 

1 12 

116 

113 

117 

1 14 

118 

115 

1 19 

1 16 

120 

117 

121 

1 18 

122 

1 19 

123 

120 

12A 

121 

125 

122 

126 

123 

127 

124 

128 

125 

129 

126 

130 

127 

131 

128 

132 

129 

133 

200 

134 


135 


136 


137 


138 


139 

201 

1 AO 


14 1 


142 


143 


144 


145 


146 


147 


148 


149 


> 



C: PROGRAM CALCULATES AZ AND EL FROM GHA AND DEC * * * ANCI AUX 
DIMENSION X( 2* 6) 

STRING OBSERC 8 ) * DATEC 15) 

F0RMATC32H ENTER LATITUDE IN DEG EG 056*3/) 

FORMATS F5 • 1 > 

FORMAT( 33H ENTER LONGITUDE IN DEG EG 256*3/) 

FORMATC 3 1H ENTER INTERVAL JUMP IN MINUTES/) 

FORMATC I 2) 

FORMATC 1 1H ENTER DATE/) 

FORMATC 5A4) 

FORMATC 16H ENTER OOOOZ GHA/) 

FORMATC 16H ENTER OOOOZ DEC/) 

FORMATC 16H ENTER 0600Z GHA/) 

FORMATC 16H ENTER 0600Z DEC/) 

FORMATC 16H ENTER 1200Z GHA/) 

FORMATC 16H ENTER 1200Z DEC/) 

FORMATC 16H ENTER 1800Z GHA/) 

FORMATC 16H ENTER 1800Z DEC/) 

FORMATC 16H ENTER 2400Z GHA/) 

FORMATC 16H ENTER 2400Z DEC/) 

FORMATC 15H ENTER OBSERVER/) 

FORMATC 1 0A3) 

FORMATC 16H AZ-EL DATA FOR * A6) 

FORMATC AH GMT* 5X* 3HGHA* 6X* 3HDEC* 7X* 2HAZ* 7X* 2HEL) 

FORMATC 25H INPUT FOR GHA IS INVALID/) 

FORMATC 37H FURTHER CALCS DESIRED THIS OBSERVER?) 

FORMATC 23H TYPE 1 IF YES* 0 IF NO/) 

FORMATC I 1) 

FORMATC I 5* AC AX* F5* D) 

FORMATC 39H' FURTHER CALCS DESIRED FOR DI FF OBSVRS?) 

FORMATC 6H LONG* 3X* F5* 1* 5X* 3HLAT* 3X*F5* 1/) 

FORMATC 5H /) 

FORMATC A15) 

WRI TEC 1* 100) 

ACCEPT BLAT 
WRI TEC 1* 102) 

ACCEPT BLONG 
WRI TEC 1*117) 

ACCEPT OBSER 
WRI TEC 1* 103) 

ACCEPT INTER 
WRI TEC 1* 105) 

ACCEPT DATE 
WRI TEC 1* 107) 

ACCEPT XC 1* 1 ) 

WRI TEC 1* 108) 

ACCEPT XC 2* 1 ) 

WRI TEC 1* 109) 

ACCEPT XC 1* 2) 

WRI TEC 1* 1 10) 
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PROGRAM 4 (cont'd) 


150 


ACCEPT XC2*2) 

151 


WRITEC 1*111) 

152 


ACCEPT XC 1* 3) 

153 


WRITEC 1*112) 

154 


ACCEPT XC 2* 3) 

155 


WRITEC 1* 113) 

156 


ACCEPT XC 1*4) 

157 


WRITEC 1* 114) 

158 


ACCEPT XC 2* 4) 

159 


WRITEC 1*115) 

160 


ACCEPT XC 1* 5) 

161 


WRITEC 1*116) 

162 


ACCEPT XC 2* 5 ) 

163 


ITIME = 0 

164 


HOUR = 0 

165 


JTIME = 0 

166 


M = 1 

167 


JB = 1 

168 


JC = 1 

169 


JD = 1 

170 


JE = 1 

171 


HTEST = 600 

172 


PI = 3.1415926 

173 


CONVT = PI/180.0 

174 


CONVTR = 1 • O/CONVT 

175 


CLAT = BLAT * CONVT 

176 


CSINL = SI NC CLAT) 

177 


CCOSL = COSCCLAT) 

178 


CTANL = CSI NL/CCOSL 

179 


R 1 = BLONG + 180.0 

180 


R2 = BLONG - 180.0 

18 1 


IFCR1 .GE. 360*0. 0R.R2.LT. 0*0 

182 


GO TO 203 

183 

202 

IFCR1 - 360.0)205*204*204 

184 

204 

R 1 = R 1 - 360.0 

185 

205 

I FC R2) 206* 203* 203 

186 

206 

R2 = R2 + 360.0 

187 


GO TO 225 

188 


C: ABOVE TESTS AND SETS R1 , 

189 

203 

JE = 1 

190 

224 

ITIME = ITIME + INTER 

191 


I FC I TIME - 60)210*208*208 

192 

208 

ITIME = ITIME - 60 

193 


HOUR = HOUR + 100 

194 


I FC HOUR - HTEST) 210*209*2 09 

195 

209 

HTEST = HTEST + 600 

196 


M = M+l 

197 


JC = 1 

198 


JD = 1 

199 

210 

JTIME = HOUR + ITIME 


TO 202 


AND R2 TO UNDER 360 DEG 
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PROGRAM 4 (cont'd) 


200 IFCJTIME - 2400) 225* 225# 300 

201 C: INCREMENTS TIME AND SETS TO HRS AND MINS* TESTS FOR 6 HR 

202 C: M INCREMENTS MATRIX ELEMENTS* 300 TAKES TO END OF DAY 

203 225 I F(X( 1*M+ 1 ) - XC 1 * M ) ) 2 1 1 * 30 1* 2 1 2 

204 211 XIN = X(1*M) - 360.0 

205 C: 301 TAKES OUT FOR INVALID DATA 

206 GO TO 213 

207 212 XIN = XC 1*M) 

208 C: TEST FOR 0 TIME PERIODS AND SET INCREMENT TO 0 

209 213 GO T0(220*221)*JC 

210 220 GO TO (214*215*216*217* 217)*M 

211 214 IFCJTIME - 0000)221*218*215 

212 215 IFCJTIME - 0600)221*218*216 

213 216 IFCJTIME - 1200)221*218*217 

214 217 IFCJTIME - 18 00)221*218*219 

215 219 IFCJTIME - 2400)221*218*300 

216 218 XINCG = 0.0 

217 XINCD = 0.0 

218 JC = JC + 1 

219 GO TO 223 

220 221 GO T0C233* 222)* JD 

221 233 XINCG = CXC1*M+1) - XIN) WINTER/ 360*0 

222 XINCD = C XC 2* M+ 1 ) - XC 2*M) )*I NTER/360* 0 

223 JD = JD + 1 

224 222 XC 1 * M) = XC1*M) + XINCG 

225 IFCXC l*M)-360. 0)235*236*236 

226 236 XC 1*M) = XC1*M) - 360*0 

227 235 XC 2* M ) = XC2*M) + XINCD 

228 C: INCREMENTS GHA AND DEC AND RETURNS TO TEST RANGE 

229 C: OUT OF INCREMENTS LOOP. CALC AZ AND EL 

230 223 PSI A = XC 1*M) 

231 IFCBLONG - 180.0)226*227*227 

232 226 IFCBLONG - XC 1* M) ) 228* 227* 227 

233 228 PSIA = XC1*M) - 360.0 

234 GO TO 229 

235 227 IFCCBLONG + 18 0 • 0 ) . GT . 360 . 0) GO TO 230 

236 GO TO 229 

237 230 PSIA = XC1,M) + 360-0 

238 229 PSI = BLONG - PSIA 

239 C: ABOVE TESTS FOR LONGS UNDER 100 AND GHA NR 360 

240 CS CONVERT ANGLES TO RADIANS 

241 PSI = PSI*CONVT 

242 CXD = XC2*M)*C0NVT 

243 CSIND = SI NC CXD) 

244 CCOSD = COSCCXD) 

245 CCOSP = COSCPSI ) 

246 SEL = CCOSL*CCOSP*CCOSD + CSIND*CSINL 

247 C: TAKES OUT FOR NEGATIVE EL ANGLES 

248 IFC SEL) 203* 260* 260 

249 260 IFC SEL - 1*0)243*242*243 


- 14 - 



PROGRAM 4 (cont'd) 


250 

242 

CEL = 0.0 

251 


TEL = 9 • 999999E36 

252 


GO TO 244 

253 

243 

CEL = SQRTC1.0 - SEL**2) 

254 


TEL = SEL/CFL 

255 

244 

CAZ 1 = CSIND/C CEL+CCOSL) - CTANL+TEL 

256 


IFCCAZ1 - 1.0)246,245,246 

257 

245 

CAZ = 0.0 

258 


TAZ = 9.999999E36 

259 


GO TO 247 

260 

246 

CAZ = SQRTC1.0 - CAZ 1**2) 

261 


TAZ = CAZ/CAZ 1 

262 

247 

ACAZ = ATANCTAZ) 

263 


ASFL = AT AN (TEL) 

264 


ASEL = ASEL*CONVTR 

265 


ACAZ = ACAZ * CONVTR 

266 


GO T0C259, 258), JE 

26 7 

259 

JE = JF + 1 

268 


WRI TEC 1, 128) 

269 


C: TEST FOR ANGLES ACROSS 360 

270 

258 

I F( C BLONG- 180.0) .LT.O.O.OK. ( BLONG+ 180*0) 

•GT. 360.0) GO TO 256 

272 


GO TO 241 

273 


C: LONGITUDES 260 TO 360 

274 

256 

I F( ( BLONG- 180*0) «LT • 0 • 0) GO TO 257 

275 


IFCR2 .GT.R1 .AND.XC 1,M) .GE.R2)G0 TO 241 

276 


GO TO 240 

277 


C: LONGITUDES 000 TO 100 

278 

25 7 

IFCXC 1,M) .LF.R1 .AND.XC 1,M) • GE* C ELONG- 180. 

GO TO 241 

28 0 


GO TO 253 

28 1 

241 

I F( BLONG - X( 1,M) ) 240, 253, 253 

28 2 

240 

I F( CAZ 1)250, 250, 251 

28 3 

250 

ACAZ = 180.0 - ACAZ 

284 


GO TO 237 

28 5 

251 

ACAZ = 360.0 - ACAZ 

28 6 


GO TO 237 

28 7 

253 

IFC CAZ1) 254, 237, 237 

288 

254 

ACAZ = 180.0 + ACAZ 

289 

237 

GO TOC 231, 232), JB 

29 0 

231 

JE = JB + 1 

29 1 


WRI TEC 1, 1 19 ) OBSER 

29 2 


WRI TEC 1, 127) BLONG, BLAT 

293 


WRI TEC 1, 129) DATE 

294 


WRI TEC 1, 120) 

29 5 

232 

WRI TEC 1, 125) JTIME,XC 1,M),X( 2, M ) , ACAZ , ASEL 

29 6 


GO TO 224 

29 7 

301 

WRI TEC 1,121) 

298 

300 

WRI TEC 1, 122) 

299 


WRI TEC 1, 123) 

300 


ACCEPT K 

301 


IFCK - 1)302,201,201 

302 

302 

WRI TEC 1, 126) 

303 


WRI TEC 1, 123) 

304 


ACCEPT K 

305 


IFCK - 1)303,200,200 

306 

303 

STOP 

307 

> 


END 

1 c 
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